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Abstract. Quantum trapping potentials for ultracold gases change the landscape of 
classical properties of scattered light and matter. The atoms in a quantum many- 
body correlated phase of matter change the properties of light and vice versa. The 
properties of both light and matter can be tuned by design and depend on the 
interplay between long-range (nonlocal) interactions mediated by an optical cavity 
and short-range processes of the atoms. Moreover, the quantum properties of light get 
significantly altered by this interplay, leading the light to have nonclassical features. 
Further, these nonclassical features can be designed and optimised. 
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1. Introduction. 

Optical lattices (OL’s) offer the ultimate control of atoms trapped by them. This leads to 
the formation of correlated phases of matter [1] , this being useful for quantum simulation 
purposes j2] and quantum information processing (QIP) applications. The degree of 
precision achieved with them so far, has allowed to achieved self-consistent light-matter 
states in a Bose-Einstein condensate (BEG) inside an optical cavity [3llll|5]. Using the 
dynamical properties of the light [6] the structural Dicke phase transition was achieved 
forming a state with supersolid features |3]. However, the study of the full quantum 
regime of the system has been limited to few atoms illlDlIIIlDSlDS]. As the light matter 
coupling is strongly enhanced in a high finesse optical cavity in a preferred wavelength, 
the atoms re-emit light comparable with the lasers used in the trapping process. As a 
consequence, an effective long-range (nonlocal) interaction emerges driven by the cavity 
held. It is now experimentally possible to access the regime where light-matter coupling 
is strong enough and the cavity parameters allow to study the formation of quantum 
many-body phases with cavity decay rates of MHz na and kHz [15]. The light inside 
the cavity can be used to control the formation of many-body phases of matter even 
in a single cavity mode PEIEIIHI. This leads to several effects yet to be observed 
due to the dynamical properties of light |T9l 1201 1211 1221 123]. Moreover, it has been 
shown that multimode atomic density patterns can emerge, even their coherences can 
become structured and light-matter quantum correlations can control the formation of 
correlated phases. Thus, a plethora of novel quantum phases due to the imprinting of 
structure by design in the effective light-induced interaction occurs m- In addition 
to light-scattering [ 13123 ], homogenous quantum many-body phases can be measured 
by matter wave scattering [261 EH EEl [29] and dynamical structure factors can be 
obtained via homodyne detection [30] • Recently, density ordering has been achieved 
with classical atoms Further, multimode cavities extend the range of quantum 

phases even further [la Ea iMiiMi. Therefore, by carefully tuning system parameters 
and the spatial structure of light, one can design with plenty of freedom the quantum 
many-body phases that emerge. The quantum nature of the potential seen by the atoms 
changes the landscape of correlated quantum many-body phases beyond classical optical 
lattice setups. Very recently an optical lattice in an optical cavity has been realised [35] 
and self-oganized Mott-insultator phases have been achieved [36] 

Moreover, the interplay between short range processes, such as on-site interactions 
and tunneling, and long-range cavity induced interactions can change significantly 
the properties of the light in the system. As these processes compete to optimise 
the energy in the system, the back-action of the matter affects the light generating 
nonclassical features um- We show how such nonclassical effects of the light inside 
the cavity arise due to the emergence of structured quantum phases of matter. This 
can be traced back to the particular structure of the full light-matter state, which 
we construct beyond the limit where the light can be integrated out (adiabatically 
eliminated) [9] EHl |37] . The formulation of the explicit form of the light-matter state 
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of the system and deriving the effective matter Hamiltonian incorporating the effect 
of light at the quantum level is a difficult problem. We provide an alternative to 
those methods that allows for the construction of the effective Hamiltonian, where 
the effect of local processes (regular atomic tunneling and on site interaction) is both 
considered in the properties of the light. The series of terms that arise, compose different 
hierarchies of light-induced interactions in addition to the adiabatic limit. This leads 
to a new effective Hamiltonian where the effect of local processes, such as tunneling 
and on-site interactions, and the global structure imprinted by the light is relevant. 
We use the technique of canonical transformations constructing the set of unitary 
operators to remove the non-diagonal terms due to the light [3H1 EH] • The underlying 
symmetries broken by design by pumping light into the system modify the structure 
of both matter and light and the competition between global and local processes are 
the origin of nonclassical features. We find the effective matter Hamiltonian and 
that the full light-matter state is a superposition of squeezed coherent states. These 
depend on the emergent quantum many-body phases of matter the system supports 
and their structural properties. We demonstrate how the quantum (quadratures) and 
classical (amplitude) properties of the light encode information about the strongly 
correlated phases of matter. As a corollary of our results, we find the conditions 
to optimise quadrature light squeezing in the system and the effect of the structure 
induced to the matter. Thus, our work will foster the design of this kind of states 
and their possible application towards quantum multimode systems in the analogous 
interdisciplinary field of optomechanics HD]. Towards possible applications, there is 
an active interest in achieving large light squeezing in optomechanical systems where 
relevant achievements have already been made [TT]. Recently, using trapped ions |12] 
superpositions of squeezed states have been achieved as proposed by |13] . Additionally, 
stationary entanglement of photons and atoms in a cavity has been studied [H], seeding 
patters via the cavity field BH and quantum control projection |l5l HO], opening the 
venue for applications on QIP. Beyond the quantum properties of light and matter, 
we find the effective master equation that describes the evolution of the system. This 
enables the possibility to study the effect of measurement back-action and its direct 
interplay with local processes. Additionally, this can be used for state preparation 
using state projection via measurement back-action ^81 091 [50], while engineering of 
non-trivial correlated quantum states is possible ISDisg, and opens the possibility to 
optimise nonclassical properties of light. 

2. The system. 

The system consists of atoms trapped in an OL inside single-mode cavity with the mode 
frequency ojc and decay rate n in off-resonant scattering. The pump light has amplitude 
Hp (in units of the Rabi frequency) and frequency Up (Ac = Up — uQ. The system is 
illuminated in a plane transverse to the cavity axis (not necessarily at 90°). The cavity 
mode couples with the atoms via the effective coupling strength g 2 = g^p/{2Aa), with g 
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Figure 1. Cold atoms trapped in an optical lattice subject to a quantum potential 
created by the light inside a single-mode cavity. The unsharp potential contour 
schematically depicts quantum fluctuations of light, which induce the light-matter 
correlations on top of the classical optical lattice potential, created by external laser 
beams. The cavity can be a standing- or traveling-wave. Different colours represent 
atoms corresponding to different light-induced spatially structured atomic modes. The 
superposition of squeezed coherent states corresponding to each light-induced mode is 
depicted on the right. 


the light-matter coupling coefficient and is the detuning between the light and atomic 
resonance [HI IMl |53]. This can be described by the Hamiltonian T-L = 1^} + 'H°‘ + 7/“^, 
where Td} is the regular Bose-Hubbard (BH) Hamiltonian [5H |55] , 

Vf’ = + h.c) - -f - 1), (1) 

(hi> * * 

with to the nearest neighbour tunneling amplitude, U the on-site interaction and pi the 
chemical potential. The light is described by 'H°' = hwcofd and the light-atom interaction 

is HH: 

= g;dF^ + g2afF (2) 

with F = D + B. D = is the density coupling of light to the atoms, 

B = E<,,„ + h-c.) is due to the inter-site densities reflecting matter-held 

interference, or bonds [211 E3]. The sums go over illuminated sites Ng, and nearest 
neighbour pairs {i,j). The operators bl (hi) create (annihilate) bosonic atoms at site 
i, d) (a) photons in the cavity, while the number operator of atoms per site is given 
by hi = is the relevant contribution to the quantum potential seen by atoms 

on top of classical OL described by the BH model, where the on-site interaction U and 
hopping amplitude to are short-range local processes. The effective parameters of the 
Bose-Hubbard Hamiltonian with the cavity held can be calculated from the Wannier 


Quantum properties of light scattered from many-body phases of ultracold atoms 
functions and are given by 


5 


^0 



VbL(x))ta(x - Xj)d''a;, 


J, 




J ta(x — Xj)M*(x)Mp(x)t(;(x — Xj)d"'x, 


( 3 ) 

( 4 ) 


where are the cavity and pump mode functions and tc(x) are the Wannier 

functions. The classical optical lattice potential is given by VbL(x). The classical 
optical lattice dehning the regular Bose-Hubbard Hamiltonian is weakly dependent of 
the cavity parameters. The atoms are mainly trapped by the strong classical lattice, 
which is created inside a cavity by external laser beams. This external potential is 
insensitive to the quantum state of atoms. The light scattered into the cavity constitutes 
a quantum perturbation of the strong classical potential. This perturbation strongly 
depends on the many-body atomic state. Additional nonlinear dependence on the Bose 
Hubbard parameters leading to semi-classical effects can be incorporated as in [191 EO] . 
Further, the classical optical lattice and cavity light can be detuned from each other. 
The light is pumped from the side of the main axis of the high Q cavity, at an angle not 
necessarily at 90°. The system is depicted in FigjT| where the effect on the scattered 
light is shown and will be explained through the paper. 

Moreover, it is useful to exploit the spatial structure of light as a natural basis 
to dehne atomic modes, as the coupling coefficients Jtj can periodically repeat in 
space [211 EQI eh E3]- atoms equally coupled to light belong to the same mode, 
while the ones coupled differently belong to different modes (p. Then we have for the 
atomic operators, 

F = + (5) 

r •r' 


where the light induced “density” and “bond” 


mode operators, such that: 


Np = and 

ietp (i,j)&ip 


( 6 ) 


with corresponding to the posible values of Ji^ and Jb,^' corresponding to Jij 

where the pair {i,j) are nearest neighbours. For example, when illuminating in the 
diffraction minima JD,ip± = and JB,ip± = FJb where Jd/b are some constants. 
These encompass the different sets of values taken by the Wannier overlap integrals 
Ji,j [2ll EH EH EH- addition, it is useful to define new emergent mode structures 
Je,p corresponding to JF = Jjj — Ji^, with {i,j) nearest neighbours. In this case in 
diffraction minima, Je,ip± = ±27^^. Note that these emergent bond terms are absent for 
structureless light, that is, scattering as it happens in the diffraction maxima of light 
{Ji,i = const). 
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3. Method 

3.1. Hilbert space rotations. 

In order to describe the physics of the system, one can construct an effective matter 
Hamiltonian from H. However, difficulty arises because the operators F and the BH 
Hamiltonian do not commute in general. Eliminating the light in the adiabatic limit, 
one can construct the effective Hamiltonian by different methods 01201EZ]. However, 
the description is only accurate as long as the magnitude of the detuning Ac is very 
large compared with any other energy scale. Beyond the adiabatic limit, our method 
relies on a series of canonical transformations constructed to eliminate the non-diagonal 
terms from the light-matter Hamiltonian. In our method, the cavity decay rate k has 
been introduced phenomenologically to reproduce the limit of adiabatic elimination of 
the light held. We hnd the additional corrections due to the non-commutativity the 
light-matter interaction with the local processes of the matter part of the Hamiltonian. 
Essentially, we perform a sequence of rotations on the Hilbert space using the formula, 

“ In.R] 

n = exp(-/?)Hexp(i?) = H + V ^ (7) 

n\ 

n=l 

with [77, the n-th order commutator with respect to R. A sequence of rotation 
operators R is performed, where these are chosen to remove the non-diagonal part of 
the light held in the light-matter interaction via the commutator expansion after their 
action. The particular structures needed and the number of rotations depend on the 
underlying structure of the matter part Hamiltonian in the light-matter interaction and 
the BH model. Their interplay with the by-products of each rotation determines the 
consecutive rotation operator to be constructed. Thus, one provides an ansatz for each R 
operator and eliminates according to the by-product of the next transformation [22] • The 
result of our method is a theory that incorporates in a perturbative operator expansion 
the interplay of the non-commutative character between local processes (tunnelling and 
on-site interactions) and the long-range (nonlocal) light induced effective interactions. 

3.1.1. Adiabatic limit rotation. The first transformation recovers the limit when light 
can be adiabatically eliminated, explicitly this is: 

-Rad = cFaf — c*F^d. (8) 


Therefore, 


g-Rad ^ exp(cFa'*' — c*F^d) = D{cF) 


( 9 ) 
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is analogous to the displacement operator from quantum optics [71 EDI EU, c = 
( 72 /(Ac + in) is the cavity Purcell factor. Thus using ([^, 

H' = 

5^eff / Af A , T^T^\\ I 


= - A.a^a + V} + + FF^) + [F, F^] 

+ ['H'',F]ca' - [W‘,F*]c‘a 
+ h[M‘,qca* - [H‘’,F<]c‘i,F]ca' 

- h[H‘,F]ca*- [H\F^]c'i,F']c'i + ... 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


to remove the additional non-diagonal light terms beyond the adiabatic limit ( |IT| ) we 
need to perform additional rotations as we state below. Dynamical effects occur due to 


(12) and light squeezing originates from (13) and (14), while the dots refer to higher 


order photon processes due to the commutator expansion that will not be considered in 
what follows. 


3.1.2. Rotations due to dynamical effects. A subsequent series of rotations is employed 
to eliminate term by term the emergence of non diagonal terms in the light sector due 
to the interplay with short range processes. The sequence of rotations is, 

Pi" = (15) 

where each rotation is given by, 

Rx,k = U’^icXkd^ + c*xld), (16) 

with k G Zq , where, = Ucjn + ^o/n, such that 




(17) 

k=0 ^ ^ k=0 ^ ^ 

Additionally, we cast our results using the natural choice of dimensionless expansion 
parameters, which are dehned as: to = to/Ac and U = U/Ac, while 

^ ^ JB,(fiSk,ip and iji^ ^ ^ JE,(pJk,(p^ h even, 

ip p 

Zk = '^ JB,p,Jk,^ and ^ JE,p>Sk,p>, k odd, (18) 

V ‘f 

where we have used collective weighted “bond” operators and weighted “bond 
current” operators Jk,ip corresponding to the light induced modes (p. These are dehned 
as, 

^ (19) 

dj>e</3 

5 ^ Anl/blb^+PX) 

dj>ev3 


( 20 ) 
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with Ahij = fij — fii. Thus these collective operators are spatially modulated by the 
difference in density between nearest neighbour pairs (?, j). The above operators can be 
traced back to the modulations to the densities and currents that arise due to the fact 
that short range tunneling and/or on-site interactions do not commute with light-matter 
interaction components. The particular binomial structure of Qn and /„ arises as each 
rotation operator needed to diagonalize the Hamiltonian generates higher order operator 
polynomials terms recursively due to the commutator expansion. These operators 
correspond to the emergence of particle-hole excitations and matter self-interactions in 
analogy to the Feynman diagram expansion in momentum space. Considering t^ and 
U as expansion parameters it is enough to consider the first few terms in the expansion 
of order (0(-)) linear in to, U and their product. Thus and are the leading 
terms, however the full expansion can be used depending on the coupling strengths and 
partial re-summation of families of terms can be employed as in the standard techniques 
of many-body physics [38l 139] . 

3.1.3. Squeezing rotation. The final transformation to obtain the effective matter 
Hamiltonian of the system is achieved by, 

= exp[(ro= - pr-)l2] = S(t) (21) 

which is the analogous squeezing operator from quantum optics [71 EDI El]- The 
squeezing amplitude operator is given by, 

I = -'fuY (22) 

<p' ^ 

where it is useful to define “density fluctuations” operators lS.N^p and “bond current 
fluctuations” A 

^4 = E (23) 

These operators of quantum fluctuations of matter arise due to higher order light 
processes entangling the matter and the light. The transformation is used to eliminated 
non-diagonal terms due to higher order photon processes. We have restricted the 
expansion to two photon processes. Therefore this contains the first nontrivial 
correction. The effect of additional higher order photon processes could be included 
in principle performing additional rotations, but we will not pursue this in what follows. 

After straight forward algebraic procedure applying each transformation and 
computing the relevant commutators using Q, we get the effective atomic matter 
Hamiltonian: 

-Hgg = ... g^x.OgfiAd-^g-RAdg-Rx.o... g-«x,oog-ife _ (^2A) 

where, Ra = —AcCffa, which is the light part after the rotation to the pump frame of 
reference [9]. 
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4. Effective Hamitonian and Full light-matter state. 

4-1- Effective Hamiltonian. 

The effective atomic Hamiltonian after the rotations is 


^eff — ^ad + + Tdy, 


(25) 


with, 


= «'■ + t ^ ptp) 


(26) 


the result in the adiabatic limit of light [21] with |k/Ac| -C 1, |H/Ac| 1, |to/Ac| -C 1 
and (/eflf = Ac| 5 ' 2 p/(A^ + nf) = Ac|cp. Here we have neglected the non-abelian shift in 
0 , but in general it can be non-zero depending on the light mode functions for complex 
values i.e. travelling wave configurations where via either pump or cavity modes light in- 
between sites can be focused. As it has been shown [211 El] , this leads to the formation 
of structures of density and bond modes that can be nearly independent from each 
other. It is possible to generate in a single mode cavity spatial multimode structures of 
R density modes [SO] EH] and 2R bond modes [21] by carefully choosing how the light 
is pumped into the system [53]. In the adiabatic limit (26), the structure of matter 


is controlled by the interplay between the BH processes, regular atomic tunneling and 
on-site interaction, and the light induced interaction proportional to The ground 


state of (26) will be achieved whenever atoms scatter light maximally for g^.^ < 0 or 


minimally > 0 m. As we will show below, the additional terms are related to 
light squeezing and arises due to the dynamical corrections from the light induced 
processes and their interplay with the short-range BH processes due to the structure 
imprinted on the matter. These will modify the landscape of quantum phases the system 
can access, as well as, the properties of light beyond being a superposition of structured 
coherent states Using the light induced mode decomposition, then we can write, 

ptp + ppt = +d 


B ,B Q Q 




(27) 


with + C.C.), where we have used the light induced “density” and 

“bond” mode operators. 

The additional contributions in "Heff are the first order corrections in tg and U 
originated by the non-commutative nature between local processes and the global 
structure introduced due to the light induced modes. The terms due to light induced 
dynamics are 


Rx ~ desip 






5'eff 


n=0 


( 28 ) 
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These dynamical terms can have a significant effect on the effective Hamiltonian 
renormalizing the effective light indnced interaction terms. Additionally, they can 
aid the formation of strnctnred gronnd-states dne to their dependency on the atom 
nnmber difference between nearest neighbonr sites. In general, these terms tend to 
induce strnctnre in the atomic density as \U\ increases dne to light-matter qnantnm 
correlations ges > 0 and dne to semiclassical effects for ges < 0 m- In principle, going 
beyond the pertnrbation character of the expansion conld be handled via renomalization. 
Close to a strnctnral phase transition, where (Ahjj) ~ 0, for example, from a strnctnred 
ground-state (with DW order) to a homogenous ground-state (a normal superffuid), the 
leading behaviour shows an instability for \U\ < 1 as. 


'^U^^'ixiXn + XnXn. 


n=0 


xl xo + Xox l 
1-U^ 


(29) 


The formation of this instability means a structure ground state can be an energetically 
favourable solution depending on the coupling constants strengths and the competition 
with other processes, from the adiabatic limit and the BH model. This provides an 
amazing potential for manipulation with the purpose of quantum simulation, as one 
can select the inhibition or enhancement of the interplay with local processes. One can 
design this using the structure constants Jjj and ges, thus controlling the light induced 
mode formation and changing the onsite interactions via Feshbach resonances or even 
via the classical optical lattice potential. The expansion could be further manipulated 
by the use of diagrammatic tools, including well know partial resumations techniques 
(i.e. Feynman diagrams) but we will not pursue this here, as we are interested in the 
regime where fo and U are perturbation parameters. The terms due to the interplay with 
short-range processes and contain the effect of higher order correlated processes, 
as they contain in addition to 2-point correlations, n-point correlations with n ^ 4. 

Beyond the adiabatic limit and dynamically generated processes, the additional 
terms in Hefi modify the energy due to light squeezing are: 

= 9esio ^ + AJ^). (30) 

ip' tp 

AN^ and AJ^ are strongly smeared out in the limit where the effect of the light is 
classical. This occurs because the atoms maximise light scattering to reach the ground- 
state of the effective Hamiltonian. However, these are relevant for the case where the 
strong classical signal is suppressed and the effect of light-matter quantum correlations 
is significant PI- In general, they have a suppression effect upon fluctuations for 
U > 0 while they promote an instability for H < 0, as gesU = \c\'^U. In addition, 
they introduce a direct modification to tunneling amplitudes that are controlled by the 
pattern of emergent bond modes via the different possible couplings Je,^- Importantly, 
as we will show their origin has a non-trivial effect on the properties of light. 
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4 . 2 . Full Light-Matter state. 

As shown previously, the effective Hamiltonian Hefr is diagonal in the light sector and 
first order quantum corrections have been included. It follows that the full solution to 
the light-matter state can be written as: 

I*) = (31) 

where the subscript “a” (“6”) corresponds to the light (matter) part; r^(f) = 
exp(—zHefft), and = r^|(pq)b. The light components are squeezed coherent 

states |a,^)a = Zl(a)S'(a)|0)a, with the squeezing operator S{f) = exp[(^*a^ — faf‘^)/2] 
and the displacement operator D{a) = exp(Q;a'*' — a*a) [H [ 6 OI |6l]. The ground 
state of the effective Hamiltonian is IT);, = \Fq)b- The light amplitudes due to 
the projection of the matter structure are ot^f\ipq)}y = cF\Lpq)b, a'^^\Lpq)b = cf}\ipq)i,, 
with fj = U^Xn- The weights due to the dynamical character of the light 

are = exp(|cp X]r=o with Xn,q>^Fq)b = Xn\Fq)b- In addition, the 

squeezing parameter amplitudes corresponding to the projection onto the matter sector 
are = i\Fq)b with, the squeezing amplitude operator Therefore, the structure 

of the strongly correlated matter gets imprinted in the quantum properties of light via 
the squeezing parameter projections \g^q)h- This generates a non-trivial superposition 
of squeezed coherent states entangled with the strongly correlated matter. In the 
above, we have neglected next-nearest neighbour and higher processes, considered all 
the first order corrections O(to) and terms order while constraining to see the 

effect of two photon processes. However, this is not a limitation in our method, since 
additional n-neighbour processes or additional n-photon processes can be incorporated 
straightforward if relevant. 

5. Quantum properties of the scattered light. 

5.1. Photon number. 

The number of photons can be written as. 

Nr 

{of a) = |c<^j^(sinh(r^j2 + \a^f) 

q=l 

= {smh{\i\r) + ^^{&G + GG^) ( 32 ) 

where G = F + f) and the weights corresponding to the matter 

component projections (the probabilities), are the coherent state 

components and the corresponding squeezing parameters depending on the 

projections on the matter states, and the number of light induced components. 
Therefore the quantum properties of the matter are accessible at the level of the photon 
number. Moreover the light-amplitude is sensitive to the particular structure that 



A(X-)-(dB) 
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Figure 2. (a) Quadrature components for light-induced two mode density coupling. 
Quadrature components are centred around a± (red and blue), the quadrature 
is centred around zero (purple). The quadrature widths are given by AXq « 
e’'(l -I- 2[dp - ^ sinh(r) cosh(r)/(l -|- ^)]e“^’')^/^/2 and AX- = e“’'(l -I- 2^[dp -|- 

sinh(r) cosh(r)/(l -|- )^/^/2, with r the squeezing parameter, dp = (Tp(A(h+)^ -I- 

A(n_)^) and Cp = (^effAcXgJ|,)/(A^ -|- k^), with |k/Ac| ^ 1. The system 
exhibits squeezing in the quadrature (j) = 7r/2. Each coherent state component is 
squeezed at the same angle, (b) Quadrature squeezing at ^ = tt/2 in dB, |k/Ac| = 
0.25(blue),O.l(green), 0.05 (purple), dashed lines correspond to squeezing projecting 
to one component, solid lines to the full state, (c) Optimal Squeezing as a function 
of the cavity decay rate, optimal squeezing is achieved whenever r = ln(|A/K|)/2 
solid line. Dashes correspond to projecting to one component with two light induced 
modes, dots corresponds to projecting to a single component with 4 light induced 
modes. Parameters are: dp = 1 (b) and (c). 

emerges due to the correlated phases of matter [21]. In the limit of large detunning 
I Ac| 3> {to, U} the above reduces to, 

+ (33) 


which is equivalent to the adiabatic limit |24j . 
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5.2. Light Quadratures and Sgueezing. 

The squeezing in the quadratures can be written as, 

. Nr 

^2sinh(r^j2) 

<?=! 

^ Nr 

- - ^ |c^J^ cos(6' - 20) cosh(r(^) sinh(r^) 


X]|c^,r(e 


\2 

Vq/ 


KqWe 


with = (e-Aa + eAat)/2, A{X^f = (X|) - {X^)\ and 9 = arg(2mAe + - A^). 

This can be rewritten as, 

A(W)^ = i + t(smh(|t|)^> 

- ^ cos(9 - 2(#.){cosh(|5|) smh(|5|)> 

+ ^A(G)^+^A(Gt)^ 


+ !T((GtG + GGt>-2(Gt>{G)) 


In particular when a^p = cG^ and d* = c*Gp so that G is Hermitian, then: 


A(Tt=) + ^(sinli(lll)t 

- ^ “s(9 - 2(#.){cosh(|5|) smh(|5|)> 


5fefr(AcCos(0) - fi;sin(0))^ 


A(G)^ 


A^iAl + nQ 

For two light induced modes with density coupling 7 ^ 0) in mean-field 

approximation, the above reduces to: 

- 2 ^ sinh(r) cosh(r) 


A(Ao)^ ^ 


1 -I- -A 

^ A2 


Alw; 


+ ^(A(n+)" + A(fl_)" + 0 (f„)) 
6“^’’ f sinh(r) cosh(r) 

~ : 1“ WG7 ^ 


1 -I- ^ 
J- A2 


+ “^ ( ^ j (A(h+)^ + A(h_)^ -h O(to)), 


(38) 
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where r = |(.^)|, and cXp = (ges^cNsJo )/The additional terms to regular 
squeezing are due to atomic fluctuations in each light induced component. We have 
used the identities, 


g±2r ^ 

= sinh(r)^ ± sinh(r) cosh(r) + 


(39) 


COS (^arg {{y + if) + = ± “ l) ’ (^0) 

and that in mean-fleld approximation A(iV+ — N_f ^ Ns{A{n_f + A{n^f) [24] . 
Beyond mean-fleld approximation, as the system goes in the normal SF state {\to/U\ 3> 
0) additional corrections due to additional coherent amplitudes that depend on their 
atomic fluctuations will increase the super-poissonian character of the light quadratures. 
The relation between quadrature components is shown in Figj^. In particular, the 
quadratures at (j) = 11/2 can be squeezed several dB in the case where the cavity 
detunning is of the order of the recoil energy as in [15] with \k/A c\ ~ 0.1 — 0.25 
and |Ac| ~ 1 — lOOF'/j. Choosing ~ Er/Ns and cTp ~ 1 with p = 3/2 in the SF 
state where fluctuations are maximal, one can easily achieve r ~ 1.5 and about lOdB of 
squeezing improving the ratio |fi:/Ac| 0.05 by changing the detunning or improving 
the cavity, see Figig^. Currently, squeezing with microwave fields has reached lOdB EH- 
Minimising ( [M| , optimal squeezing is found when r = ln(|Ac/K|)/2. Optimal squeezing 
at 0 = 7r/2 is given by. 
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+ f (A) (A(ftn= + + 0((„)) , (41) 


Therefore, optimal squeezing is limited strongly by how small is the ratio \k/A c\. Note 
that each coherent state component is squeezed stronger by a factor of 1/i?, squeezing in 
each component ~ 50% more with respect to the total state, see dashed lines in Figj^b) 
and (c). Thus the projection to a single component improves optimal squeezing. This 
projection is even natural, being a consequence of spontaneous symmetry breaking in 
the system without optical lattice [63] . In the case of R density modes one has in general 
for |k/Ac| < 1, 
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(42) 


with hq corresponding number operator of the light-induced density mode component 
per site. In general, projecting to a single component of R modes produces an 
enhancement factor on squeezing of 1/i?. For example, projecting to a single component 
with four light induced modes gives a enhancement factor of 75%, see Figj^c) dotted 
line. Therefore, by incrementing the number of light-induced modes one can optimise 
squeezing in a single component even though the cavity decay ratio \k/A f\ is not that 
small. The general structure of the light-matter properties for arbitrary number of 
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bond and density modes is rather involved as it contains information regarding the 
correlated phases of matter that emerge. It is instrnctive to see the effect on the 
sqneezing parameter r for some cases, as we will show in what follows. 

6. Squeezing parameters and emergent strnctnred phases. 

6.1. Homogenous light scattering. 

When atoms scatter light homogeneonsly {JB,ip = Jb, Jd,^ = Jd, JE,ip = 0), local 
density imbalance is snppressed. As it has been shown [la EH [IZ!, in the adiabatic 
limit SF and MI scatter light differently depending on the properties of the qnantnm 
many-body state. As other energy scales become relevant, the additional terms amonnt 
to renormalization of the indnced interaction in the matter wave coherences, the 
terms in BB\ snch that = gesU'^B'^/Q — U"^) for |f7| < 1. This allows to 
enhance the effect dne to the matter wave coherences via the on-site interaction in 
the effective Hamiltonian. The sqneezing parameter using mean-held approximation is 
r = 1(^)1 = 2z\geEU\J%Ns\'ilJ*‘^ip‘i) + + {fif) - 2n‘f - ni\/Al, with atom 

number per site Ui = (hi) and the SF order parameter = {hf). Therefore, for a 
deep classical OL {.Jb = 0) there will be no squeezing due to r. As the classical OL 
becomes shallower {Jb 7 ^ 0), the squeezing parameter is maximal for MI and smoothly 
decreases as we reach the SF state, see Fig. [^. This is correlated with the fact that 
light scattering while illuminating in between density maxima (at the bonds) is maximal 
in the MI while decreasing as the SF grows [5H] . 

As the number of light induced modes in the matter increases, the induced 
structures play a substantial role on light squeezing. For 2 light induced modes, such 
that their amplitudes alternate sign every other site {Jd,ip = '^Jd or Jb,,^ = ±Jb and 
Je,ip = 2 J£)), we hnd that the matter induces structure to the squeezing parameter. 
As it has been shown [21] besides from SF, MI the system supports gapped superfluid 
states, dimer phases, supersolid (SS) and density waves (DW). 

6.2. Diagonal coupling, illuminating at lattice sites. 

Without bond ordering {Jb,^^ = 0, Jo,^ = ^zJb), the squeezing parameter is different 
for SS, DW, SF and MI phases. In mean held theory, the squeezing parameter is 
proportional to the product of SF order parameters in each light induced component, 
r = 2z\ges\toJDNs{'iljf'if_ + c.c.)/Af where V'lt correspond to each light induced mode 
component. Thus, for an insulating state (DW or MI) r = 0 while for a SF = '0_ 
and in the SS state -0+ 7 ^ fj-. Indeed, as the onsite interaction increases i.e. for 
half integer hllings when atoms scatter light maximally to reach the ground state 
of Hes (fi'eff < 0), light will be squeezed maximally in SF, while as SS emerges, it 
will diminish until reaching the DW state where no squeezing is possible, see Fig. 
S). The total coherent state amplitude is 7 ^ 0 when DW order is present, while 
di = 0 in the SF or MI. Moreover, when atoms scatter light minimally to reach 
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Figure 3. Squeezing parameter across different transitions, (a) Scattering with a 
single light induced mode {Jd,lp = 0, Jb,o = JB,geS > 0), as the system goes from MI 
to SF the squeezing parameter diminishes while the change becomes more visible as 
density increases p = l.O(solid), 2.0(dashed), 3.0(dotted)). (b) Two component system 
for maximum light scattering {Jd,± = ^Jd, Jb,ip = 0,5eff < 0) at p = 3/2 filling. The 
system goes from DW to SS and to SF as zto/U increases, dashed line is the difference 
in density Ap = |p_|_ — p_| and solid line is the squeezing parameter r. (b) Inset, MI 
to SF transition {Jd,± = ^Jd,Jb,ip = 0,peff < 0) at p = 1. (c) Squeezing parameter 
r for the supersolid dimer (SSD) to SF transition {Jd,ip = 0,Jb,± = < 0) at 

p = 1.0. For zto/U < 1 the system is in SSD while for zto/U > 1 is SF. Parameters:(a) 
Peff = U/Ns, Jb = 0.05, f = rAl/{2z\gef(U\JlN,) (b) peff = -0.5t7/iV„ Jo = 1.0, 
f = rAl/{2z\gesto\JliNs); (c) peff = -25U/Ns, Jb = 0.1, f = 2rAl/{z\geeU\J'/^Ns). 
Ns = 100, z = 6. 

the ground-state {pes > 0), the squeezing parameter is different for SF and gapped 
SF states. The squeezing parameter for a quantum superposition (QS) state [2ll HE] 
is tqs = 4:z\gef[\toJf)Ns{m + l)(nj — m)(l + m — nf)/Af for incommensurate fillings 
m < Hi < m + 1 with m positive integer. Thus, for a gapped SF r < uqs while for a 
normal SF state r > rqs- 

6.3. Emergent bond order. 

In addition, dynamical terms can induce bond ordering due to the emergent coupling 
JE,ip as U increases. Emergent bond ordering due to density coupling occurs because 
products of weighted bond and bond current operators modify the effective Hamiltonian 
via These terms arise because on-site interaction and tunneling do not commute 
in general with the light-induced long-range interaction. The new terms that appear 
in the effective Hamiltonian favour density imbalance as \U\ increases and modify the 
coupling of matter wave coherences with it. Explicitly, we have to order 0{tyj‘^), 

^ 2g,^Jltl[{l - lJ^){ClJo, + h.c.) + U\b\B^ + h.c.)] (43) 

The current operators Jq = are structureless but the weighted bond operators 

= E(i — hi)(blbj + bjbi) induce a staggered field between bonds 

as density varies between every other site and the difference between atom populations 
can alternate sign. Essentially, the density variation acts as an additional dynamical 
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diffraction element that affects the interference of the matter waves in between density 
maxima. The matter waves in other to compensate the staggered field and optimize the 
energy in the effective Hamiltonian acquire a phase pattern between adjacent sites. This 
translates in the formation of dimer states. The difference in phase of the matter waves 
A0 7 ^ 0. Thus, for < 0, bond ordering will occur and dimer physics [21] will emerge 
even in a deep optical lattice. As a consequence, 4 bond light induced modes will form 
leading to a superposition of 4 light-matter correlated squeezed coherent states. The 
squeezing parameter can be cast as r = z\ges\toJD^s[4>i+4>3+i4>3+4>4) cos(A0)]/A^ with 
4>q = ^0 = arg(-^ 2 ) ” = arg(' 04 ) — arg(-^i) and fjq the order parameter 

of each effective induced mode. Therefore, the squeezing parameter of light inherits the 
structure due to bond-ordering even in a deep optical lattice. Therefore, the interplay 
between short-range processes and the long-range cavity induced interaction leads to 
the emergence of physics absent in the classical optical lattice and the adiabatic light 
limit in this configuration. This implies a new alternative for the design of Hamiltonians 
containing dimer physics in analogy with spin-liquid Hamiltonians |59j . 

6.4- Off-diagonal coupling, illuminating in between lattice sites. 

In the case with only off-diagonal light-matter coupling {Jd,ip = 0, = ±Jb, 

Je,ip = 0 ) for maximal light scattering (^fes < 0 ), the squeezing parameter is different 
between SF, superfluid dimer (SFD), supersolid dimer (SSD) and SS states. The 
effective interaction strength in the adiabatic limit gets renormalized by the term 
-|- i?fi?))/2 to order 0{U‘^). The squeezing parameter in mean- 
field approximation and using the typical statistical properties of the states [H 2 | can 
be estimated as r zlgesUlJfNglcpl + 03 + (02 + 01) cos(2A0) — n^ns|/(2A^), where 
the populations for each dimer are ha and hb- In contrast to diagonal coupling 
7 ^ 0), the normal SF state r is minimal (r —)■ 0), as 0i = 02 = 03 = 04 = |0P, 
nA = nB = 2no ~ 2|0p and A0 = 0. In SFD, bond ordering occurs, thus 0i = 03, 
02 = 04, nA = ub = 2no, A0 7 ^ 0 with r 7 ^ 0. Typically, dimer states have 
7 r /2 < A0 < 71, thus cos(2A0) < 0 depending on the parameters chosen for the system. 
For SSD, bond ordering and density modulation occurs, then 0i 7 ^ 03, 02 7 ^ 04, ua 7 ^ ub, 
A0 7 ^ 0. Thus, r is maximal as uaUb < 4no, see Fig. i:. Therefore, as bond ordering 
occurs and DW order emerges r is different from zero. For minimal light scattering 
> 0) one has direct information regarding SS order solely due light-matter quantum 
correlations. The squeezing parameter has considerably simpler structure with respect 
to dimer phases since, r z\gef[U\JfNs{n^ — n_)^/A^ with n± the number of atoms in 
each light induced mode per site. The squeezing parameter tracks directly the emergence 
of DW order. Therefore, for SS r 7 ^ 0 and for homogenous SF r = 0, while the coherent 
state amplitudes are a± = 0 but the number of photons is {d'^d) 7 ^ 0. 
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7. Effective Master Equation. 

Beyond the |k/Ac| ^ 1 limit, we use the methods of quantum optics [501 El] and we 
hnd the effective master equation for the system as 

^ = — ^[HeSyP] + (2G^pG' + [&G,p]^) (44) 

P ^ifiq)a\'pq)b bipilai^ipiy 

‘Pq,Vl 

where p is the density matrix, [•,•]+ is the anti-commutator and Pq^i are the matter 
coefficients (probabilities) that can describe either a pure or mixed state. The second 
term in the master equation is the effective Liouvillian which includes dissipation. 
Measurement back-action beyond the |fi:/Ac| 1, |17/Ac| -C 1 and |to/Ac| -C 1 limits 
can be devised by using G as the effective jump operators for quantum trajectories. The 
operator G is related to the coherent state amplitudes of light via = G|4/). In 

the above the Markov approximation is implied in the limit when T = 0 [601 El]- This 
allows to consider the effect of measurement back-action, the role of local processes and 
their interplay due to light-induced non-local interactions simultaneously. This opens a 
new venue for exploration regarding the design of global structured dissipation channels 
and measurement induced projection and state design of non-trivial quantum correlated 
states [SniEDEl], as well as, control [IHl EE] and the transition to classicality [El] . As it 
has been shown, this can greatly enhance and optimise the desired quantum properties 
of light by design. 

8. Conclusions. 

We have shown that quantum optical lattices offer a new tool to engineer a generalised 
class of states that are a non trivial superposition of structured squeezed coherent 
states of light entangled with matter. These states are entangled with the matter at 
the fundamental level due to the structure of quantum many-body matter states. We 
have demonstrated that breaking symmetries by design one can induce structure to the 
parameters that control the nonclassical features of light. This has been shown to be 
accessible via quantities such as, the photon number and the quadratures of light. We 
have shown how the quantum properties of light contain the information of matter-field 
coherences, density patterns of matter and light-matter quantum correlations. Thus 
the properties of strongly correlated phases of matter get imprinted on the quantum 
properties of light. Moreover, we have found that the interplay between induced long- 
range processes and ordinary short-range atomic processes lead to the modihcation of the 
effective Hamiltonian of the system. We have obtained that the effect of local processes 
can be used to generate delocalized dimer phases due to the dynamical properties 
of light even in deep optical lattices. Moreover, one can optimise these nonclassical 
features depending on cavity parameters and the structure imprinted to the matter 
that gets transferred to the light-induced mode structure. Additionally, the non-trivial 
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light-matter correlated states that arise can be used to design dissipation channels 
via the effective master equation and measurement back-action. A pathway to study 
the behaviour we describe is to combine several recent experimental breakthroughs: 
detection of light scattered from ultracold atoms in OL was performed, but without a 
cavity [B51 EB] and BEC was trapped in a cavity, but without a lattice [31 in E]- As 
our treatment of the system is based on off-resonant scattering, this is not sensitive to 
a detailed atomic level structure. Therefore, our treatment applies to analogous arrays 
of natural or artihcial quantum objects such as: spins, fermions, molecules (including 
biological ones) [67|, ions [68], atoms in multiple cavities [69|, semiconductor UDl or 
superconducting qubits HH. 
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